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f**^ ■ Field quantization in unstable optical systems is treated by expanding 

£ — ' the vector potential in terms of non-Hermitean (Fox-Li) modes. We define 

f*"^ , non-Hermitean modes and their adjoints in both the cavity and external 

regions and make use of the important bi-orthogonality relationships that 
exist within each mode set. We employ a standard canonical quantiza- 
tion procedure involving the introduction of generalised coordinates and 
Q,' momenta for the electromagnetic (EM) field. Three dimensional systems 

are treated, making use of the paraxial and monochromaticity approxi- 
mations for the cavity non-Hermitean modes. We show that the quantum 
C^ ' EM field is equivalent to a set of quantum harmonic oscillators (QHO), 

associated with either the cavity or the external region non-Hermitean 
modes, and thus confirming the validity of the photon model in unsta- 
ble optical systems. Unlike in the conventional (Hermitean mode) case, 
the annihilation and creation operators we define for each QHO are not 
^^ ' Hermitean adjoints. It is shown that the quantum Hamiltonian for the 

J"j , EM field is the sum of non-commuting cavity and external region contri- 

butions, each of which can be expressed as a sum of independent QHO 
Hamiltonians for each non-Hermitean mode, except that the external field 
Hamiltonian also includes a coupling term responsible for external non- 
Hermitean mode photon exchange processes. The non-commutativity of 
certain cavity and external region annihilation and creation operators is 
associated with cavity energy gain and loss processes, and may be de- 
scribed in terms of surface integrals involving cavity and external region 
non-Hermitean mode functions on the cavity-external region boundary. 
Using the essential states approach and the rotating wave approxima- 
tion, our results are applied to the spontaneous decay of a two-level atom 
inside an unstable cavity. We find that atomic transitions leading to cav- 
ity non-Hermitean mode photon absorption are associated with a different 
coupling constant to that for transitions leading to photon emission, a fea- 
ture consequent on the use of non-Hermitean mode functions. We show 



that under certain conditions the spontaneous decay rate is enhanced by 
the Petermann factor. 



1 Introduction 

The photon model of the quantum electromagnetic (EM) field has recently been 
developed further for the areas of cavity QED and quantum optics in dielectric 
media. Recent work has demonstrated its validity in microscopic theories Q, 
where EM field quantization is the same as in free space, and in macroscopic 
theories, where the classical optical device or material medium is treated as 
a spatially dependent permittivity. Canonical quantization carried out by ex- 
panding the vector potential in terms of the true (or universe) modes for the 
system (obtained from a Hclmholtz equation) shows that the EM field is again 
equivalent to uncoupled quantum harmonic oscillators (QHO), one for each true 
mode ||, ||, [|, U . If idealised or quasi modes (as obtained from an approximate 
optical system permittivity) are used instead of true modes, coupled QHO arc 
obtained, one for each quasi mode |3|, g] and photon exchange processes be- 
tween quasi modes occur. Phenomena such as beam splitter effects |7j, light 
energy loss from cavities jq], reflection and refraction can be explained via 
such processes using quasi mode theory. Both localised quasi mode (as in cavity 
or external regions for the standard cavity QED model M) and non-localised 
quasi mode (as in the beam splitter 0) cases occur. Quasi mode theory has 
been extended (lO) to cases || where two sets of quasi modes are needed. 

The photon model obtained by quantization based on normal mode (free 
space, true or quasi mode) expansions involves mode functions determined from 
Hcrmitean eigenvalue equations and which satisfy standard orthogonality con- 
ditions (power orthogonality). However, for unstable optical resonators or gain- 
guided lasers, the natural EM field modes to use (Fox-Li modes) are eigenfunc- 
tions of a non-Hermitian operator. Biorthogonality conditions apply for these 
non-Hermitean modes (NHM) and the related set of Hermitean adjoint mode 
functions (HAM). The appropriateness of the photon model Jo] has been ques- 
tioned for such systems, where phenomena such as linewidth enhancement in 
gain-guided lasers occur. This linewidth enhancement was interpreted |12| as 
excess spontaneous emission into the laser mode and described by the Peter- 
mann factor, and several semiclassical studies (for example |13| [TJ, Q) have 
successfully accounted for the experimentally observed Peterman factors using 
the normalisation integrals for the HAM functions (the original NHM are nor- 
malised to unity) . A variety of differing quantum treatments aimed at explain- 
ing the linewidth enhancement have been published. Some have the objective 
of trying to account for the success of the semiclassical theories in terms of a 
fully quantum treatment of the field, others present other explanations of the 
excess noise without using the NHM functions. In p6[ the enhanced linewidth 
is attributed as a propagation effect arising from longitudinally inhomogeneous 
noise gain, rather than an excess of local spontaneous emission, the microscopic 
rate of spontaneous emission into a given non-power-orthogonal cavity mode 



not being enhanced by the Petermann factor. Reference |l7f accounts for the 
excess noise factor via either spontaneous emission noise or amplification of vac- 
uum modes leaking into the cavity, depending on the choice made for operator 
ordering. This treatment did not involve the use of NHM and single atom decay 
was found to be enhanced by the cavity Q factor, not the Petermann factor. 
In Q a simple few mode model was presented, and excess noise attributed 
to loss-induced coupling between the cavity eigenmodes, with no use of NHM 
being made. In reference |19| on the other hand, field quantization based on 
the NHM and the HAM is described, and the spontaneous emission rate of an 
atom in the cavity found to be enhanced by the Petermann factor. The work 
in |2C[| is a laser theory involving a reservoir of excited atoms. Starting from a 
true mode description, quasimodes and their adjoints satisfying non-Hermitean 
eigenvalue equations are introduced, and results approximating to those for the 
semiclassical theories are obtained for the noise amplification. Finally, in refer- 
ence [£1[ field quantization based on NHM and the HAM obtained by solving the 
Helmholtz equation in both the cavity and external regions with non-Hermitean 
boundary conditions is presented. As in |17J the single atom decay was found 
to be enhanced by the cavity Q factor, not the Petermann factor. 

The present paper is a further quantum treatment involving a standard 
canonical quantization proceedure based on expanding the vector potential in 
the unstable cavity region via the non-Hermitean (Fox-Li) modes and their 
Hcrmitean adjoint modes. It is aimed at trying to account for the success of 
the semiclassical theories in terms of a fully quantum treatment of the field. 



Our approach is similar to those of Q, 21 , although here we employ the 
canonical quantization method and also treat three dimensional systems, using 
the paraxial J22] and monochromaticity approximations J19[ for the cavity NHM. 
As in |l9f, plfl , we consider all field modes, rather than a small number as in 
the simplified model treated in 111]. Our NHM and HAM are also obtained 
from the properties of the unstable optical system, rather than from treating 
atom- field interactions, as in [J20| . As in |21|], the external region is treated 
and the field described via further NHM and HAM. Similarly to [|||, 0, we 
find that the field is equivalent to a set of QHO's, associated now with non- 
Hermitean modes rather than true modes and thus confirming the validity of 
the photon model for the case of unstable optical systems. Similarly to fll9f , 
fell , the annihilation, creation operator pairs for each QHO are not Hcrmitean 
adjoints. However, our results differ significantly in detail from those in |LS|, fl2lf . 
Our canonical quantization procedure involving both right and left travelling 
modes leads to a doubling of the annihilation, creation operator pairs compared 
to Efl], pl[ . We show that the total number of true mode QHO's equals the 
total number of QHO's associated with either with the cavity or the external 
region NHM, one oscillator for each annihilation, creation operator pair. As in 
[Ell, the final form of our quantum Hamiltonian for the EM field is the sum 
of non-commuting cavity and external region contributions, but is simpler as 
there are no off-diagonal terms. To a good approximation, both the cavity and 
external region Hamiltonians can be expressed as the sum of independent QHO 
Hamiltonians for each NHM, but the external region Hamiltonian also includes 



a coupling term responsible for external NHM photon exchange processes. The 
two independent QHO Hamiltonians are alternative choices for an unperturbed 
Hamiltonian. Analogous to ]19[ for the cavity region Hamiltonian, we obtain 
left and right eigenstates for each of these unperturbed Hamiltonians and the 
energy is given by the usual QHO result. As in plf, certain cavity NHM and 
external region NHM annihilation, creation operators do not commute, leading 
to cavity energy gain and loss processes, though again the details differ. We 
find a simple description of the resultant cavity-external region light coupling 
in terms of surface integrals involving products of cavity and external NHM 
functions. Atom-field interactions are treated in the electric dipole and rotating 
wave approximations. We show that atomic transitions leading to cavity NHM 
photon absorption are associated with a different coupling constant to that for 
atomic transitions leading to photon emission. This feature is directly related 
to the treatment being based on NHM functions and leads to enhanced emission 
rates. Using the essential states approach, we consider spontaneous decay of a 
two level atom located in the cavity. The coupling term in the field Hamiltonian 
is neglected, assuming that atomic decay into the cavity is much faster than 
cavity decay. Coupled equations are obtained for atom-field states amplitudes 
involving the excited atom and no photons or the ground state atom with one 
photon in a cavity NHM. Markovian decay occurs under certain conditions, the 
decay rate being enhanced by the Petermann factor p2] in special cases when 
a single NHM dominates. 

The plan of this paper is as follows. In Section @ quantization of the EM field 
based on NHM is carried out and approximate energy eigenstates determined. 
In Section atom-field interactions are considered and a simple treatment of 
spontaneous emission from a two level atom into an unstable optical cavity is 
presented. The results are summarised in Section 0. The Appendix contains 
the detailed calculation of certain commutation rules. 



2 Theory of canonical quantization via non-Hermitean 
modes 

2.1 Non-Hermitean mode functions 

The case of interest is that of unstable optical cavities, for which a typical 



example is shown in figure 1. Using the paraxial approximation |22|, 23, |24j , 
the cavity NHM functions U„(R) and HAM functions V„(R) describing right 
propagating fields associated with light of angular frequency uj n can be written: 

U„(R) = Afa n .exp(ik n z).u n (s,z) (1) 

V n (R) = Afa n .exp(ik n z).v n (s,z), (2) 

with the polarization vectors chosen as a n = (x or y) and the wave number is 
k n = LJ n /c (> 0). A/" is a normalization factor. The fields u n (s,z) and v n (s,z) 
are slowly varying functions of z describing the transverse s = (x, y) dependence 



of the NHM functions. These fields satisfy the right propagating paraxial wave 
equation: 

(Vy + 2ik n d/dz) u„(s,z) = (V§. + 2ik„d/dz) v n (s,z) = 0, (3) 

where Vt = x d/dx + y d/dy is the transverse component of V. In addition, 
u n (s,z) and v n (s,z) satisfy eigenvalue equations of the form: 

£ u n (s,z) = 7„ u„(s,z) (4) 

£} v n (a,z) = 7* v n (a,z), (5) 

where £ is a two dimensional non-Hermitean round-trip propagation operator 
describing right travelling light, £^ is the related Hermitean adjoint operator 
related to the transpose operator £ T describing left travelling light via £' — 
(£ T )*, and r y n are the complex eigenvalues which are independent of z. The 
operators £, £ T and £^ are linear integral operators that transform a field in 
the z plane into a new field in the same plane, and are defined in terms of the 
propagation kernel G+(s,z; sq,z) such that for any arbitary transverse function 
f{s,z) we have: 

g(s,z) = £ f{s,z) = d 2 s G + {s,z:s ,z)f(s ,z) (6) 

h(s,z) = £ T f(s,z) = / d 2 s G + {s 0l z;s,z) f{s ,z) (7) 

i(a,z) = 2 f f(s,z) = / d 2 s G* + (s ,z;s,z) f{s ,z). (8) 



Thus £ maps f(s,z) onto g(s,z), £ T maps f(s,z) onto h(s,z) and £' maps 
f(s,z) onto i(s,z). The propagation kernel G+(s,z;sq 7 z) is obtained from the 
Fresnel approximation to the Huygens-Fresnel principle, details are given in 
references E3, G4fl. For the situation depicted in figure 1, G+ would be con- 
structed allowing for right travelling propagation from the z plane to the right 
end mirror, left travelling propagation to the left end mirror and then right 
travelling propagation to the z plane. The fields £ f(s,z) and £^ f(s,z) also 
satisfy the same paraxial wave equations as u n (s,z) and v n (s,z). In our notation 
n will specify & n , k n (the polarization and longitudinal wave number) and the 
transverse mode index 9 n for the distinct transverse NHM functions in which 
the quantities & n , k n are the same. It will be convenient to quantize the wave 
numbers, with k n — N n n/ I , where N n is an integer and / is the cavity length. 
Both U„, V„ describe right propagating light fields, whereas 

W„ = a n .exp(ik n z).w n (s,z) = (V„)*, (9) 

describes left propagating light fields, where w n (s,z) = v n (s,z)* and satisfies the 
eigenvalue equation £ T w n (s,z) — j n w n (s,z), . Thus we see that u n (s,z) and 
v n {s,z) are special solutions of the right propagating paraxial wave equation that 
are self-reproducing under the propagation operators £ and £' respectively. 



Based on the well known two dimensional orthogonality results |l3f for the 
u n and v n such as J d 2 s u n (s,z)*v m (s 1 z) = S nm , we can obtain important three 
dimensional orthogonality integrals in the form: 

/d 3 RU„(R)*.V m (R) = S nm (10) 

Jc 

f d 3 RU„(R)*.U m (R) = Cum (11) 

Jc 

/<i 3 RV„(R)*.V m (R) = D nm , (12) 

Jc 

the first being the biorthogonality condition, the others giving the overlap in- 
tegrals for the NHM and the HAM in terms of the Hermitean, positive definite 
transformation matrices C and D. The integrals are over the cavity region. It 
can be shown that C nm and D m „ are zero unless the angular frequencies uj n , u> m 
and the polarizations & n , & m are equal, and C D = D C = E. By convention 
we choose C nn = 1, and then from Schwarz' inequality it follows that D„„ ^ 1. 
The Petermann factor K n is given by D„ n in the semiclassical theories (see 
[ [Hj[ [14] , ]l5[). The NHM functions satisfy a completeness relation: 

£(CC(R)Kf (R7 + C(R)*^(R')) = S afi 6(R - R) (13) 

n 

for R, R' inside the cavity, and where U° (R) , V™ (R) are the a components 
(a, (3 — x,y) of the vector fields U n (R),V n (R). Integrals similar to those in 



equations ( JlOj , 11, ^2|) but without the complex conjugation are ignored due to 
the product of the two fast-varying exp(ifc„z) factors approximately averaging 
to zero. The corresponding completeness relationship in reference [fill does not 
include the second term. There are also interrelationships between the NHM 
and HAM functions in terms of the transformation matrices: 



U„(R) = J2 C "- V ™( R ) V «( R ) = X! D ™» U -( R )- ( 14 ) 



2.2 Vector potential and generalised coordinates 

The NHM functions and true mode functions represent monochromatic fields. 
For more general time dependent fields we can write the vector potential A(R) 
as the sum of a right travelling light field A^(R) and a left travelling field 
Ar(R) for R inside or outside the cavity: 

A(R)=A H (R) + A R (R)*, (15) 

and then expand A#(R) as a linear combination of suitable monochromatic 
fields. One choice is to expand A# in terms of true mode functions Afc(R) 



in the right half space R. Corresponding true modes in the left half space are 
denoted Afc*(R) = (A fc (R))*. These modes all have the same angular frequency 
Uk ■ Alternatively, a similar expansion can be made (but only within the cavity) 
in terms of either the NHM functions U„(R) or the HAM functions V„(R). The 
expansion coefficients (qk or Q n or R n for the three choices) act as generalised 
coordinates specifying the vector potential in the appropriate region. In the left 
half space qk* = {qk)* ■ Thus: 



A fl (R) = J2q k A k (K) (16) 

k 

A fl (R) = ^Q n U„(R)=]Ti?„V„(R), (17) 

n n 

where R lies within the cavity in equations (|17|). The Q n , R n depend on the qk, 
as the latter specify the field everywhere via equation (fig). We have: 

R R 

Qn = 2^ r„fc q k R n = 2_^, ^ nk ?*" (18) 



k k 



where 



T nk = f d 3 RV„(R)*-A fe (R) (19) 

Jc 

A nk = /rf 3 RU„(R)*-A fc (R). (20) 

Jc 

The matrices T,A have properties T f A = A f r = E and (A t A) nTO = C nm , 
(r'r) nm = D nm . Interrelationships between the Q n and the R n can easily be 
found in terms of C, D. These are: 



(°c n / *-*nm -^m **>n / *-'nm Wm V J 

771 771 

Being determined from solutions of the paraxial wave equation, the cavity 
NHM and HAM functions U„ and V„ actually represent the lowest order terms 
in a more general treatment starting from paraxial fields. Following Ref. p2| , 
the spatial part of the vector potential in the Coulomb gauge 

V-A = (22) 

for harmonic solutions of the wave equation (and which represent right travelling 
waves), can be expressed as 

A^" } (R) = exp( i fc n z).[F T (R) +zF z (R)], (23) 



where Ft and z Fz are transverse and longitudinal components. Ft and Fz can 
be expressed as power series in the small quantity / = Wq/Iu, where wq is the 
lateral size for the beam and Ik = k (wo) 2 is the Fraunhofer diffraction length. 
Ft,Fz involve even, odd powers of / respectively. The zeroth order term in 

Ft satisfies the paraxial wave equation p2|, thereby clarifying the status of 

— . . (k ) 

this equation as determining the lowest order contribution to An* . The cavity 

NHM and HAM functions U„ and V„ are particular zeroth order solutions for 

Ft- However, to satisfy the Coulomb gauge condition even in zeroth order, the 

first order term from Fz would need to be included. This would replace U„ 

and V„ by: 

U« = (U n +i/k n z{V T -U n }) (24) 

V« = (V n +^nz{V T -V n }), (25) 

resulting in the gauge condition being satisfied correct to the second order in /. 
In a higher order treatment we would replace U ra and V n by these expressions, 
which now would have a small first order correction in the z direction. However, 
it is apparent that the small correction terms are perpendicular to the zeroth 
order terms U n or V„. This results in expressions to terms in the Lagrangian 
and Hamiltonian based on NHM expansions only changing in second order, 
and such corrections can be neglected in our results. Although the Coulomb 
gauge condition plays a central role in the field quantization presented here, the 
correction terms in the NHM functions need not concern us further, since our 
formal canonical quantization procedure is based on the true mode functions 
[|j. These are constructed to satisfy the Coulomb gauge condition for the free 
space situation that applies both inside and outside the cavity. 

Expansions of the right travelling field A^(R) may also be made in the ex- 
ternal region in terms of NHM Ujf(R) and their HAM Vjj-(R) as determined 
from other non-Hcrmitean operators, and the expansion coefficients Qk,Rk 
act as generalised coordinates for the field in the external region. We do not 
need to specify the actual method for determining these mode functions other 
than noting that the paraxial approximation would not apply. Generalisations 
to three dimensions of the approach in plfl might be applied. However, the 
external NHM functions should satisfy the Coulomb gauge condition and the 
biorthogonality conditions, and be chosen so that the magnetic energy contribu- 
tion to the Lagrangian is in diagonal form (see equation (EO)). Orthogonality, 
completeness and interrelationships analogous to equations (HI O, 111, Il3, pTl) 
occur for the \5k (R) and *Vk (R) ■ Equations analogous to equations (Il7 , Qq, 
^l[ ) apply for the Qk,Rk- The matrices C, D, T, A are replaced by G, H, A, $, 
where: 

fd 3 RU K (R)*.V L (R) = S KL (26) 

J E 

f d 3 RUif(R)*-Ui(R) = G KL (27) 

J E 

/d 3 RV^(R)*-V L (R) = H KL (28) 

J E 



and 

®Kk 



[ d 3 RV K (R)*-A k (R) (29) 

JE 

f d 3 RV K (H)*-A k (R). (30) 

JE 



The integrals are over the external region. The matrices G,H, A,$ satisfy 
relations analogous to those for C, D, T, A. 

The vector potential in general satifies the wave equation. Since the Q n 
or R n determine the vector potential inside the cavity, they also determine A 
and its derivatives on the boundary with the external region. Hence they also 
determine the held in the external region, which is specified by the Qk or Rk ■ 
However, the field in both regions is also specified by the q k . Hence the sets Q n 
or Rn, Qk or Rk are equivalent to the q k as generalised coordinates specifying 
the vector potential, and these coordinates are interconvertable. The inverse 
relationships are: 

% = E M uKl Qn = £ Nj^ Q K = Yl M ^ T nl #» = £ N M A h #* 
l,n l,K 1,71 l,K 

(31) 
where the matrices M, N are defined by cavity and external region integrals: 

M kl = /d 3 RA t (R)*.A,(R) (32) 

Jc 



>c 



/ d 3 RA fc (R)* -A^R). (33) 

JE 



IE 

Clearly Mj,; + N k i = S k i .The relationship in equations (|3^) illustrates the point 
that the number of true modes k in the half space must equal the number of 
cavity NHM designated n or the number of external NHM designated K . Hence 
the total number of true modes in the full space would be equal to twice the 
number of cavity or external NHM. Equivalently, the total number of true modes 
Afc in the full space would be equal to the number of cavity right travelling NHM 
U„ plus left travelling NHM W„ = V*, or the number of right travelling XJk 
and left travelling V^ external NHM. This result is important in comparing the 
total number of quantum harmonic oscillators for the field described via true 
modes with that when it is described by NHM. 

2.3 Lagrangian 

The Lagrangian is defined by: 

L = ¥°J d3R (A(R)2 " C2[V X A(R)]2) - (34) 



The true mode expansion gives L in terms of the generalised coordinates q k and 
velocities q k (and their complex congugates) in the form: 

R 

L = £v^{qkq* k -ulq k q* k ). (35) 

k 

The Lagrangian can also be written as: 

L = L C + L E , (36) 

the sum of contributions Lc, Le from the cavity and external regions. Using the 
NHM expansion the cavity contribution Lc can be obtained in terms of either 
the generalised coordinates Q n and velocities Q n or the alternative quantities 
R n and R n (and their complex congugates). The external contribution Le 
involves Qk, Qk or Rk, Rk- The expressions for these Lagrangians are: 

Lc = e J2(QnK-uj 2 n Q n R* n ) = eoJ2Mn-"l R nQ*n) ( 37 ) 

n n 

Le — e o / , {QkR*k ~ u k QkR*k) = £ o / t (RkQk ~ u k RkQ*k)- (38) 
k k 

For the cavity term the magnetic field integral has been evaluated approxi- 
mately with derivatives of the slowly varying factors being neglected. The ex- 
ternal NHM and HAM have been chosen so that the matrix Q, determining the 
magnetic energy contribution, which is given by: 



/d 3 R[VxU K (R)]*-[VxV L (R)] (39) 

Je 



n KL = c 2 

I E 

u 2 k 5kl, (40) 



is in diagonal form with eigenvalues ui K . For the external NHM a wave number 
&&- can be defined via ojk = cfcx- 

2.4 Generalised momenta 

The true mode generalised momenta are defined in terms of the Lagrangian and 
related to the generalised velocities as: 

Pk = dL/dq* k (41) 

Pk = £o9fe- (42) 

In the left half space p k * — ijpk)* ■ We can also write p k as the sum of cavity 
and external terms as: 

Pk = Pk+Pk (43) 

p c k = dLc/dql pZ=dL E /dql (44) 
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For the cavity NHM it is convenient to introduce generalised momenta defined 
in terms of the cavity term in the Lagrangian via: 



P, 



dL c /d(Q n y s n = dL c /d(R n y 



(45) 



rather than expressions involving derivatives of the total Lagrangian L. We are 
able to make this somewhat arbitary choice for P n and S n because our formal 
quantization process is based on the true mode quantities pk and qu, and so 
only the pk need to be defined in the canonical way. The generalised momenta 
P n and S n are given by: 



Pn — £0 / y ^iim Qm — ^oRn 
m 

Sn = £0 / y t)„ m JX m = £oVr 



(46) 
(47) 



We note that the generalised momentum P n associated with the generalised 
coordinate Q n is given in terms of the generalised velocity R n associated with 
the alternative coordinate R n . The reverse applies as well, and these features 
are a direct consequence of the two mode functions U„(R) and V„(R) not 
being the same for a non-Hermitean mapping operator £ , since u n ^= v n . 
For the external region similar considerations give for the generalised momenta 
Pk = coRk and Sk = £oQk- Together the sets of NHM generalised coordinates 
and momenta Q n , Pn or R n , S n or Qk,Pk or Rk, Sk are equivalent to the true 
mode generalised coordinates and momenta qk,Pk- The cavity P n ,S n depend 
on the pk- We have: 



Pn — 2^ Kik Pk S n — 2_^ Tnfe Pk 



(48) 



Interrelationships between the P n and the S n can easily be found in terms of 
C, D. These are: 



^n / ^nm ^m ^n — / *-*nm *rr 



(49) 



Equations analogous to equations (Q |9|) also apply for the Pk, Sk- The overall 
inverse relationships are: 



Pk 



E M *X^ 



£ N kl^Kl fit = £ M M A nl Sn = £ N U*h ^ 
l,K l,n l,K 

(50) 
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2.5 Conjugate momentum field 

The conjugate momentum field n(R) is defined as the derivative of the La- 
grangian density with respect to A and is given by: 

II(R)=e A(R). (51) 

For all positions R it can be written as the sum of a right travelling light field 
Il^(R) and a left travelling field Ilij(R)*. The generalised momenta acts as 
expansion coefficients for expressing the conjugate momentum field in terms of 
the mode functions. In terms of the true mode and NHM expansions involving 
the generalised momenta we find that: 

(52) 

(53) 

^P„V„(R) (54) 



n(R) = 


= n fl (R) + n fl (R; 




R 


k(R) = 


-- ^ Pfe A fc (R) 

k 


[fl(R) = 


-- Vs„U n (R) = 



with R inside the cavity in equations (54). Analogous expressions apply for 



the external region. Comparing the similar expression for the vector potential 
(equations (|l7j)) we note the unexpected feature that the momentum S n is now 
the expansion coefficient for the NHM function U n (R), whilst P n is that for 
the HAM function V„(R). This role reversal for the generalised momenta is a 
consequence of introducing the bi-orthogonal NHM functions. 

2.6 Hamiltonian 

The Hamiltonian H is defined via the true mode expression and obtained in 
terms of true mode generalised coordinates and momenta: 

R 

H = X>*3fc+P*&]-£ (55) 

k 

H = Yl ( pkP *k + e o w fc*?fc)/ 2e o- (56) 

fc 

The Hamiltonian is that for independent harmonic oscillators, one for each true 
mode in each half space. To obtain the Hamiltonian in terms of NHM we first 
write the Hamitonian as the sum of cavity and external contributions He, He 

H = H C + H E (57) 

given by 



He = J2\pU*k+Pk*4k]-L c (58) 

k 
R 

He = J2^q*k+Pk*qk]-L E . (59) 
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We then show that 

Pk = Y,( A nkQn + K k Rn)/2eo (60) 

n 

Pk = Y,^KkQK + A* Kk R K )/2e . (61) 

K 

Equations ( pL8[ ) etc. then enable us to eliminate the k sums involving <jk,<ik, 
so that He, He now only contain NHM generalised velocities Q n ,Rn,QK,Rn- 
Eliminating all the NHM generalised velocities in favour of the NHM gener- 
alised momenta via P n = eoR n , S n = eoQ n etc. finally gives the contributions 
to the Hamiltonian He, He from the cavity and external regions. The cavity 
contribution He only depends on Q n , P n or R n , S n . and the external contribu- 
tion He only depends on Qk,Pk or Rk,Sk- Expressed in terms of just one 
set of generalised coordinates and momenta (the Q, P or the R, S choice), both 
He and He are Hamiltonians for coupled harmonic oscillators, two oscillators 
for each n and K and there will be double sums over n and K . However if a 
mixture of both choices are used, the cavity and external region contributions 
to the Hamiltonian can be expressed as single sums: 

H C = ^2({Pn*S n + Sn*Pn}+e 2 U 2 n [Q n Rn* + RnQn*])/2€0 (62) 

n 

H E = Y.^ PK * SK + SK * P ^ +e l Uj2 K\QKRK*+R K QK*])/2e . 
K 

(63) 

Using the transformation relations in equations (21, 49), the sums over n of the 
two terms in each of the square brackets in equation ( p2|) can be shown to be 
equal, however we have chosen to leave the Hamiltonian in a more symmetrical 
form regarding the two sets of generalised position and momentum coordinates. 

2.7 Canonical quantization and commutation rules 

Canonical quantization involves replacing the generalised coordinates and mo- 
menta by quantum operators satisfying standard commutation rules. The re- 
placements are: 

Ik -> qk, q* k -*■ <?fc f = Qk* ,Pk -> Pk, P* k -* Pfe f = Pk* (64) 

for true modes and 

^cn y ^cn, \}°cn) * ^cn , *n ^ ^n, (,-*- n) y *n , 

Rn * Rn, \Rn) * Rn , "Jn * &n, \&n) * "n (65 j 

for the cavity NHM. The standard commutators that are non zero are: 

[9fc,P/ t ] = [qk^Pi] = ih5u (66) 
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from which the commutators for the cavity NHM operators can be obtained via 
the relationshps in equations (Qq, W8) and the properties for the T, A matrices. 
For the cavity NHM operators the non zero results are: 



[Qni + m \ — [Qn i Pm\ — ["ru <~>m J — [Wn i ^m\ — ifi nm (u7 j 



and 



[Rnj Pm \ = *" ^nm |^n > •* mj = *" ^mn- ("") 

Similar replacements are made for the external NHM generalised coordinates 
and momenta, and commutation rules analogous to those for the Q n , P n , R ni S n 
are obtained for the Qk , Pk i Rk , Sk ■ The commutators involving both cavity 
and external NHM operators will be discussed below. The Hamiltonian H be- 
comes the operator H and the vector potential and conjugate momentum fields 
A(R),n(R) become field operators A(R),n(R). The formal expressions for 
these operators in terms of the generalised coordinate and momentum opera- 
tors are exactly the same as in the classical expressions, so need not be repeated 
here. 

The commutation rules for the field operators A(R), fl(R) can be obtained 
via their true mode expressions and are given by: 

[A Q (R),fl /3 (R / )] =ihd a0 6{Il-Tl / ). (70) 

Those for the electric displacement field operator D(R) = — n(R) and the mag- 
netic field operator B(R) = V x A(R) can then be obtained as: 



[B Q (R),D (R 7 )] =ihJ2 £ e tl 3 1 d/dX 1 5(R-R / ) (71) 

7 

The commutators involving both cavity and external NHM operators are less 
straightforward to determine than those just involving one type at a time. One 
suitable approach involves using the commutation rules obtained for the electric 
displacement and magnetic field operators. The details of this calculation are of 
some interest, and are set out in the Appendix. The non zero results can be ex- 
pressed in terms of surface integrals over the boundary S between the cavity and 
external regions , and are given below as I nK , J n K ,K nK , L nK , J nK , J nK , K, nK ,£ n K- 
Apart from others obtained by taking the Hermitean adjoint, the non zero com- 
mutators are for the case of cavity generalised coordinates and external gener- 
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aliscd momenta are: 



[Qn,P 1 K ] = (h/k n ) [ d 2 sV* n -U K =ihL nK (72) 

Js 

[RnJlc] = (h/k n ) f d 2 s\J* n -V K =ihK nK (73) 

Js 

[Q n ,S f K ] = (tl/k n ) f d 2 sV; i -V K = ihJ nK (74) 

Js 

[Rn,P]i\ = {h/kn) ( d 2 s\J n -V K = ihi nK (75) 

Js 

The non zero commutators for the case of cavity generalised momenta and 
external generalised coordinates are: 

[P n ,& K ] = {ih/k 2 K ) f d 2 sz-XJ* n x(Vx V K ) = -ihJC nK (76) 
Js 

[S n ,PJ K ] = (ih/k 2 K ) f d 2 sz-V* n x(Vx\J K ) = -ih£ nK (77) 
Js 

[P n ,$ K ] = {ih/k 2 K ) f d 2 s z ■ XJ* n x (V x V K ) = -ihl nK (78) 
Js 

[Sn,Q f K ] = {ih/kjc) f d 2 sz-V; i x(VxV K ) = -ihJ nK (79) 
Js 

These commutators are much smaller than those involving only cavity quantities 
or only external quantities. 

2.8 Annihilation and creation operators 

In view of the link between the various classical Hamiltonian terms and har- 
monic oscillators, it is appropriate to replace the generalised coordinate and 
momentum operators by annihilation and creation operators . For the true 
modes we define: 

dk = Afe^fe + i/J. k pk a\ = \ k ql - i^kP\ 
d k * = Xk% + i/^kPl «L = ^kd_k - i^kPk, (80) 

where Xk = \J €quj k / Iti and fi k = ^/ 'l / '2heouj k . These have non zero commuta- 
tion rules [dk, d\] = Ski and the Hamiltonian operator can be written as the well 
known form 

H = Y / (ald k + l/2).ha;k, (81) 

k 

corresponding to independent quantum harmonic oscillators, one for each true 
mode in both half spaces. For the NHM a somewhat unusual choice is made 



15 



owing to the juxtaposition of the Q, R operators and P, S operators in the 
Hamiltonian contributions He, He- We combine the Q with the S and the R 
with the P along with the Hermitean adjoints to define annihilation and creation 
operators. 

For the cavity NHM it is possible to define eight operators of the annihilation 
and creation operator type: 



A n 


- AnC^ji + l/l n b n 


A n - 


- A n K n l\X n r n 


Bn -- 


= X n Rl + ifi n P^ 


B* = 

n 


- Antyn ^Mti n 


A n ; 


KiQn ~ l ^n^n 




= An-Rn + i/J- n P r 


Bt - 


- A n ri n l[A n ± n 


J - J n 


= KQl + in n &, 



(82) 

which are expressed as linear combinations of the eight original generalised posi- 
tion and momentum operators and their adjoints. With these choices we see that 
the annihilation, creation operator pairs are now (A n ,A$), (B n ,B#), (Af^ , A^) 
and (B*^,Bl) rather than the expected (An,Al), (£„,!?+), {A*\A*) and 
{B%< ,-Bjf). The non zero commutation rules are: 



[i„, A*] = [B n ,B*] = [A*\Al] = [B#\Bl] = S nm (83) 



and 



[AnJU = D„ m [B#\B#]=T> mn (84) 

[Bn,,B m \ = C mn [A n , A m \ = C nm . (85) 

Some commutators are zero because of the selection rule w„ = w m for non zero 
C nm and D nm . Similar commutation rules occur in reference [fill (their a n 
corresponds to our A n , their b' n to our A%), except that here we have twice as 
many annihilation, creation operator pairs since no operators corresponding to 
our B n and Bf occur. The precise nature of the eight operators can be seen by 
expressing each of them as a linear combination of the true mode annihilation 
and creation operators. Each turns out to actually involve both types of true 
mode operator, for example A n is a linear combination of the a& and the a kse 
with k in the right half space. However, the expansion coefficient for one type of 
operator involves the factor f-(Lu n ,uJk) = (u) n — u^) / (2^u) n u)k), the other the 
factor / + (aj n ,(jjfc) = (w n + u>k)/(2y/ui n uik)- Assuming that all the field modes 
of interest have approximately the same frequency (monochromaticity approx- 



imation, cf. 1 19 ) terms involve the first factor can be ignored and the second 



factor approximated as unity. We then find the approximate relationships: 



16 



A n = J]r nfc a fe A*^J2K k al 





A' 


A 


B n = 


R 

- 2^ A nfc «fc* 


R 

B#=^n k al 

k 


M" = 


- 2^ A n fc at 


At^ T nA 

k 


4 #t = 


- z2 r * nk Sfc * 


R 

k 



(86) 



We thus can identify A n and B n as annihilation operators associated with right 
and left travelling modes respectively, and with A^ and B# as their correspond- 
ing creation operators. A similar identification of A^ and B^ as annihilation 
operators, A n and B n the corresponding creation operators also applies. Using 
the selection rule lo„ — w m for non zero C nm and D„ m , the eight operators for 
the cavity NHM can be interrelated as follows: 



m m 

A\ = Y^T> mn A* B*i=J2VmnB m . (87) 



These relationships are used to simplify the cavity term in the Hamiltonian He ■ 
A similar approach can be used to construct annihilation and creation op- 
erators for the external NHM. For the external NHM we write: 



*k = 


= X K Qk + i^K^K 


A% = X K R]c - WkPr 


B K -- 


= Aa'-R^ + i^K P K 


B% = X k Qk - i^K^K 


4 - 


= X K Qk~ iHK&K 


AP = X k Rk + i^ K p K 


B K - 


= X K Rl< - i^K P K 


Bp = Xk& k + in K & K 



The commutation rules are analogous to those for the cavity NHM operators, 
except that some further commutators are non zero because the selection rule 
luk = ^l for non zero Grl and Ukl is only approximate, and the expected re- 
sults like [Ak, A' l ] = Ukl are multiplied by the factor J + (ujk, ^l) that is only 
approximately unity. Approximate relations between the external NHM anni- 
hilation and creation operators and the true mode operators that are analogous 
to equations ( |36|) and interrelations between the external NHM operators that 
are analogous to equations d87]) ca n be obtained. These depend respectively, 
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on the monochromaticity assumption and the approximate selection rules on 
GkltHkl- Operators such as Ak and Bk are approximate linear combina- 
tions of true mode annihilation operators, with A K and B K their corresponding 
creation operators. 

2.9 Field operators 

The field operators A(R),II(R) can be now expressed in terms of annihilation 
and creation operators. In terms of true modes we have: 



A(R) = J^ ^h/2e io k {{a k + a[j A fe (R)+[4 + a k *} A fc (R)*) (89) 

k 

= ^ V h /^o^k(a k A fc (R)+4 Afc(R)*) (90) 

k 

and 

« i , 

fl(R) = Y, -V^o^k/2([a k - 4J A fc (R)-[4 - &k*} A fc (R)*) (91) 

k 

= Y -V^Wfe / 2(a k A fe (R)-4 Afe(R)*) (92) 

k 

for all positions R. Expanding in NHM functions the field operators are: 



A(R) - J2 VK/2tou n ([An + B*} U n (R)+[i# + B n ] V„(R)*) (93) 
ft ( R ) = H JV^o^„/2([i„ - B#] U n (R)-[i* - 4.] V„(R)*), (94) 

n 

for positions R in the cavity region. The monochromaticity approximation 
is not required. Comparing equations (|9[ |9l]) with equations ( p3| , |94|) again 
suggests the identification of A n and B n as annihilation operators associated 
with right and left travelling modes respectively, and also indicates that we 
have twice as many annihilation operators as in reference J19[ . Their expressions 
corresponding to equations (93, 9J) do not include the B# and B n contributions. 



2.10 Quantum Hamiltonian 

The cavity term in the Hamiltonian can be expressed in terms of the cavity 
NHM annihilation and creation operators. We find that: 
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H C = J2 ^{ (\\AtAn + AlAf] + I) + {\[B#B n + B-lBp] + I)} 

n 

= J2 ^«i (4 M f + ^) + (AiAf + ^)l (95) 

n 

= ^^ n {(i#i n + i)+(S#B n+ I)}. (96) 

n 

The first form is obtained just by inverting equations ( B2n and substituting into 
the operator form of equation (p2|). The derivation of the second and third 
forms for He uses equations (p7j) > which is based on the selection rule ui n = u) m 
for non zero C nm and D m „. The final expression for the cavity contribution 
Hq to the Hamiltonian turns out to be the Hamiltonian for a set of independent 
quantum harmonic oscillators (QHO), two for each cavity NHM index n - one 
corresponding to the left travelling NHM, the other to the right travelling NHM.. 
The lack of coupling between the cavity NHM QHO follows from the result that 
all A#A n and B#B n commute with each other. Because of the commutation 
rules, the quantities A#A n and B# B n act as number operators, even though 
they are not Hermitean. However, although the individual terms in equation 
( Pq) are not Hermitean, the overall sum giving He is. The existence of QHO will 
become apparent when the approximate energy eigenstates for the Hamiltonian 
H are obtained. 

Expressions for the field operators in the external region analogous to equa- 
tions (|9J| |j) apply however without using the approximate relations between 
the external NHM annihilation and creation operators and the true mode op- 
erators that are analogous to equations ( pq ) and the interrelationships between 
the external NHM operators that are analogous to equations (pTJ), and more 
importantly the external region contribution to the Hamiltonian He can be 
written as: 



H E = ]T hw K { {\[A*A K + P K Ap] + i) + (\[B*B K + B^Bp] + -)} 

(97) 



K 



If we now apply the relationships that analogous to equations (pjj), the last 
expression can be simplified to give an approximate form for He '■ 



H E = J2 ^k{ (A*A k + i) + (B*B K + i)}. (98) 

K 

The approximate external region contribution H E to the field Hamiltonian will 
be seen to be the Hamiltonian for a set of independent QHOs, two for each 
external region NHM K. Independence follows from all the A^Ak and BJ^Bk 
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factors commuting with each other. Again, the non Hermitean quantities A* Ak 
and BJ-Bk act as number operators. The full expression for the external region 
contribution to the Hamiltonian can then be written as the sum of the harmonic 
oscillator term H E and a coupling term Ve as follows: 

H E =H E + V E (99) 

where the coupling term V E is 

Ve = IJ2 ^ A k A V + B f K BP} ~ \ Yu huJ ^ A * k « + B * B *} (100) 

K K 

It should be noted that neither H E nor Ve are Hermitean, though their sum 
giving He is. It is also possible to write: 



where: 



He = (ff°)t + (V E )\ (101) 



(£°)t = Y, nu K{(Ai K Ap + l)+{Bl c Bp + l)} (102) 

K 



2' ^-*~« 2' 
(Ve) 1 - ^Y. hhJ K{A*A K + B*B K } 

K 

-\Y, *"*{ A k A P + B lcBp} (103) 

K 

The field Hamiltonian H is now given by: 

H = H c + H e = Hc + H e + Ve = H c + (H° e )t + (Ve)+. (104) 

Although there are similarities in that the quantum field Hamiltonian is the 
sum of non-commuting cavity and external region contributions, the final form 
of our Hamiltonian for the EM field differs from that in reference [Ell. Here 
there are no off-diagonal terms so the Hamiltonian is simpler. 

2.11 Non- commutation of cavity and external NHM op- 
erators 

The commutation rules for either the cavity NHM or the external region NHM 
annihilation and creation operators separately have been described. Since the 
field Hamiltonian is the sum of cavity and external region contributions, each 
of which depend on their own annihilation and creation operators, it is also 
important to consider the commutators between cavity and external region op- 
erators. These can be obtained from those for the generalised coordinates and 
momenta and hence depend on certain surface integrals, as we have seen in 
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equations (172 - [79[). The commutators are listed in tables 1 and 2. In these 
tables a n K = (1/2) y/u> n / ujk and b n K = (1/2) \J^k / w n . Another thirty-two 
commutation rules can be obtained by taking the adjoints of the results in tables 
1 and 2. 



[L,R] 


A K 


A* 


B K 


B* 

D K 


A n 





0>nK L n K 
+ b n K CnK 


— a n K L n K 

+ b n K C n K 





A* 

n 


- a nK K nK 
-b nK IC nK 








a n K K nK 

-bnK /C* K 


B n 


- a n K K nK 

+b n K K-* nK 








a n K K n K 
+bnK K-nK 


B* 

n 





0-nK L n K 
~b n K &nK 


— a n K L n K 

— bnK^-nK 






Table 1. Commutators for cavity and external NHM annihilation and cre- 
ation operators. 



[L,R] 


A K 


m 

A K 


#K 


D K 


A n 


0-nK JnK 
+ b n K JnK 








— a n K JnK 
+b n K JnK 


A* 





— O-nK I nK 
-bnK^nK 


a nK I nK 

-b nK T* lK 





B n 





— anK I nK 
+ bnK 1 n K 


a n K I nK 

+b n KlnK 





B* 


0>nK JnK 
— bnK JnK 








— a n K JnK 

— b n K JnK 



Table 2. Commutators for cavity and adjoints of external NHM annihilation 
and creation operators. 



2.12 Energy eigenstates for unperturbed Hamiltonians 

Several possible choices are available for a suitable unperturbed Hamiltonian for 
which approximate energy eigenstates can be obtained. One choice is the cav- 
ity region Hamiltonian He, another is the external region approximate Hamil- 
tonian H^ . A third choice is the external region approximate Hamiltonian 
(H^y. In the case of He the non-Hermitean quantities A^A n , Bf[B n be- 
have as a set of commuting number operators for which we can find simulta- 
neous left and right eigenstates, even though the individual number operators 
are non-Hermitean. These eigenstates are also eigenstates of He- Similarly, 
in the case of H a E the non-Hermitean quantities A^Ak, B^Bk again behave 
as a set of commuting number operators for which we can find simultaneous 
left and right eigenstates, even though the individual number operators are 
non-Hermitean. These eigenstates are also eigenstates of H^. In each case, 
the unperturbed Hamiltonian involves the annihilation, creation operator pairs 
(A n ,A#), (B n ,B# ) or (Ak,A k ),(Bk,B*) and the true mode vacuum state 
|0 >= |...0fe Ok*-- > acts approximately as a vacuum state for the NHM 
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annihilation operators A n , B n or Ax, Bk , with 

A n \ >= B n \ >= i K | >= Bjf 1 >= (105) 

(see equations (|86|)). The right simultaneous eigenstates of the number operators 
A#A n , B#B n or A*Ak, B^Bk can be constructed in the standard form as: 



\{n n },{m n } > =l[ {A f)l" iB pl" \0>, (106) 



\{n K },{m K } > =II ( K r\ " ^ l°> (107) 

where n n ,m n ,nK,mK are all positive integers. These states are also right 
eigenstates for the unperturbed Hamiltonians Hq or Hg, where (with {Np} = 
{n n }, {m n } and {N_e} = {n K }, {m K }) 



H C \{N C } > =E({N C }) |{N C }>, (108) 

i?£|{N £ } > = £({N £ }) |{N B }> (109) 

and the energy is given by: 

E({n n },{m n }) = y^ hw n (n n + m n + 1) (110) 

n 

£({%}>k}) = ^ hw K {n K + m K + l) (111) 

These are clearly quantum harmonic oscillator energy expressions. For the state 
|{"-n},{?7in} > there are n n ,m n photons associated with the cavity NHMs 
U n {R), V n (R)*, and for the state |{^i<"},{m^} > there are n,K,rriK photons 
associated with the external region NHMs Uk(R),Vk(R)* ■ As the number 
operators A#A n , B#B n or A^Ak, B^Bk are non-Hermitean, the left simul- 
taneous eigenstates of these operators are different to the right eigenstates. 
However, both the left and right eigenstates are also eigenstates of the unper- 



turbed Hamiltonian He or H^ : having the same energy as in equations (110, 



111 ). The left eigenstates are obtained in the form: 



(K},{mJ| = <0|n^= 1 ^7=, (H2) 

„ \fn n \ V m n' 

({n K },{m K }\ = <0|n (il ^T (g / ^C (H3) 
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and the eigenstates for He satisfy ortho-normality and completeness relations: 

({N C }|{N^} > =«5({N C },{N/,}) (114) 

J2 |{N C } > ({N C }| = 1 (115) 

{N C } 

in an obvious notation. The round bracket notation for the left eigenstates has 
been introduced to distinguish ({Np} | from the adjoint of the right eigenstates, 
(| {N c } >) T = < {Nc} |, whose role will be discussed below. Analogous results 
apply for the states based on the external Hamiltonian H^. Similar results have 
been obtained in reference |L9[. 

Some simple states will be specified as examples of these expressions. The 
zero photon or vacuum state is the same as the true mode vacuum. The one 
photon states are: 



I* • =i#|0> |lf>=B#|0> (116) 

*\0> K>=B#\ 



\li > =A*\0> \1*>=B*\0> 



For the cavity Hamiltonian He there are two other choices of right and left 
eigenstates, which are the adjoints of the states ({n„}, {m„}| and |{fi n }, {mn} > 
and are given by: 



\{n n },{m n }) = n (A K" (jB ^r" |0> (U7) 

n V ft V ft 

< {n n },{m n }\ =<0|II n /~l >A ■ ( 118 ) 



Both states have energies given by equation (110) and the new right eigenvalue 
equation is: 

H c |{N C }) = £({N C }) |{N C }). (119) 

The new orthogonality and completeness relationships that apply are: 

< {Nc}|{N c }) = <5({N c },{N c }) (120) 

£ |{N C }) < {N C }| = 1. (121) 

{N C } 

These results may be obtained directly by taking the adjoints of previous equa- 
tions or by using the alternative form for He in equation (pa), noting that 
A* n A%> and BjjBf* are also number operators and that A*t| >= B%*\ >= 
using equation (|86|). 

For the approximate external region Hamiltonian the situation is similar. 
Here the additional right and left eigenstates of the Hermitean adjoint {H^Y 
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are |{nif }, { m K}) and < {tik}, {m^ll, which are the adjoints of ({nx}, {wrt}| 



and \{nK},{m,K} > respectively and can be obtained from equations (113, 



107). Both states have energies given by equation (111) but the right eigenvalue 



equation now involves (H^y and is: 

(£°)t|{N E }) = £({N £ })|{N £ }). (122) 

Orthogonality and completeness relationships analogous to the cavity case apply, 
and the results are obtained similarly. 

3 Atom-field interactions and spontaneous de- 
cay 

3.1 Rotating wave approximation and coupling constants 

In the electric dipole approximation the interaction between an atom and the 
field is given by V = d • A / e Q , where d is the atomic dipole operator and the 
conjugate momentum field II is evaluated at the atomic position (see reference 
[|J). The dipole operator can be written as the sum of a component d + as- 
sociated with upward atomic transitions and a component d associated with 
downward atomic transitions, where d + = (d _ )t. The conjugate momentum 
field can be expressed as the sum of its positive and negative frequency compo- 
nents, H + ,lT~ respectively; Thus we have: 

d = d+ + cT (123) 

n = n+ + fr. (124) 

Expressions for the positive and negative frequency component for positions R 
in the cavity are: 

n+ = V - y/he w n I 2(A n U n (R)+B n V„(R)*) (125) 

Z ' I 

n 

ft- = - V V«eo«n/2(B#U„(R)+i# V„(R)*) =(ft + ) t (126) 

* — ' t 

n 

Note that although n~ = (n + )'!' overall, the same is not true on a term by term 
basis. In the rotating wave approximation (RWA) the atom-field interaction is 
given by: 

VBWA = {d+-fl + + d--fl~)/e Q . (127) 

For a two level atom (TLA) with upper state |e >, lower state \g >, the dipole 
operator components are: d = d fr± , where d is the dipole matrix element 
< e | d | g > and <r+ = |e >< g\, <7_ = \g >< e| are the upward, downward 
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transition operators. The atomic transition energy is Hujq. For the TLA located 
in the cavity at position R, the RWA atom-field interaction can be expressed 

as: 



Vrwa = <r+ E H9n A n + gl K)\ - &- E %n B* + Q* A*)], (128) 

n n 

where the coupling constants are defined by: 



g£ = -i v /w n /2fie (d-U n (R)) (129) 

gl = -Vwn/2fie (d-V„(R)*). (130) 

These coupling constants are also referred to as one-photon Rabi frequencies. 
Here we note one of the unusual features of the NHM theory - whereas an 
upward atomic transition accompanied by the absorption of a U n (R) photon 
will be associated with the coupling constant g£, a downward atomic transition 
accompanied by the emission of a U n (R) photon will be associated with the op- 
posite coupling constant g^ . The situation involving the absorption or emission 
of a V„(R)*photon is similar, with the opposite coupling constant g„ again 
being associated with emission. Thus absorption and emission involve differ- 
ent coupling constants, a peculiarity ultimately responsible for the Pctermann 
factor. 

3.2 Decay of excited atom 

The spontaneous emission problem for the TLA located in the cavity will be 
treated via a simple approach using the essential states approximation and the 
RWA. In addition, the effect of the external region term He will be neglected 
in the first instance, where we just focus on the initial emission of photons into 
the cavity NHM. The external region term He comes into play when the loss of 
cavity NHM photons and the creation of external NHM photons is considered - 
that is, cavity loss processes. As this implies, we are assuming that the cavity 
loss rate is slow compared to the spontaneous emission rate from the atom into 
the cavity. Thus we will take as our Hamiltonian: 

H t = H a + H c + Vrwa, (131) 

where H a is the atomic Hamiltonian and ignoring the zero point energy terms 
as usual. Hence our essential states will just include product states of the form: 
|e > |0 >, \g > |l£ > and \g > |1^ >. The initial state is given by: 

|*(0)>=|e>|0>, (132) 

and the time dependent Schrodinger state vector is written as: 
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|*0) > = C e (t) exp(-iwo*/ 2 ) |e > |0 > 
+ E C »( t ) ex PHK-^o/2]t)| 5 > \ln> (133) 

+ ^C n B (t)expHK- Wo /2]t)|g > |1* >, (134) 

n 

where C e (t),C^(t) and C„(t) are amplitudes in the interaction picture. Using 
the time dependent Schrodinger equation we obtain coupled equations for the 
amplitudes: 



dC e {t) 

l ; 

dt 

dt 

• jggft) 

dt 



= - E ^™ exp( ? <5„t) <# (t) + g* oq>(W B t) C r f (£)} (135) 

n 

= 5r f exp(-i<5 n t)C e (t) (136) 

= <tfexp(-tf n t)C e (t), (137) 



where the detuning is S n = loq — uj n . Note the opposite coupling constants 
in the last two equations. The one photon amplitudes C„(t),C^(i) can be 
formally eliminated to give a single integro-differential equation for the excited 
state amplitude C e (t) of the form: 



dC e {t) 



dt ,/o 



dTK{r)C e (t-T), (138) 



where the kernel K{t) is defined by: 

K{r) = Y, 2 9n9n^M^nr). (139) 

In view of the convolution integral form, the equation for C e (t) can be solved 
in the general case via Laplace transform methods in terms of the Laplace 
transform of the kernel: 

K( S )=iJ2^9n9n/(S n + is). (140) 

n 

Recalling that our generic index n specfies the angular frequency to n , the 
polarization a n and the transverse mode index 6 n , we may convert the sum over 
n in equation ( |139| ) to an integral over io n times a mode density p(u> n ), together 
with sums over a n and 8 n . Then, assuming that the coupling constants <?„ , g^ 
are slowly varying with frequency to n , we see that the kernel K(t) decreases to 
zero over a correlation time r c of order the inverse bandwidth of the coupling 
constants. Further, assuming this time scale is short compared to the decay time 
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T e of the excited state amplitude, we can then make the Markoff approximation. 
In this case: 



dC °W~{l drK(r)}C e (tl (141) 



dt jo 

showing that the probability of the atom being excited, P e (t) = |C e (i)| 2 decays 
exponentially to zero with a rate T e given by: 

r e = -2ReK(e), (142) 

where e is small. 

To evaluate r e , equation (|14|) can be used to interrelate the coupling con- 



stants in equation (14C), giving 

^ = -^D nm (^)*, (143) 



which is then substituted into equation (14C). Now the transformation matrix 
D nrra is zero unless the frequency and polarizations are the same (uj n — u> m , a„ = 
a m ), so after replacing the sums over n and m by an integral over uj n involving 
NHM density p(u> n ) and by sums over a n and the transverse mode indices 
6n,9m, we then obtain for K(e) the expression: 



K(e) = i f dio nP (co n )2j2J2T,^n)'Dn m (gir r - 

J a a l W " 



P 



(«n - Wo) 
-IT f dw n p{uj n ) 2J2Y.J2 (9n) Dnm (<?„)* %„ - W )(144) 

J ™ a a 



From the Hermiteancy properties of D nm it is not difficult to see that the double 
sum J2 e 2 e (0n ) D ram (9m)* i s rea ^ so that in general we have for the decay 
rate: 

T e = 4np(uj )Y / T,\9n(^)\ 2 K n (oJo) 

+47rp(w ) £ Yl Yl 9n^o)'Dn m (^o)9 m (^o)*, (145) 

where we have substituted the Petermann factors K n (u)o) for the diagonal ma- 
trix elements D m (uo). Note that all (diagonal and off-diagonal) transforma- 
tion matrix elements D„ m with a n = a rn are independent of a n . In general 
then, the decay rate T e will involve the cavity NHM coupling constants g£ and 
the transformation matrix elements D„ m evaluated at the atomic transition 
frequency luq . The decay rate is the sum of diagonal terms involving the Peter- 
mann factors K n (u>o) — D„ n (wo) and off-diagonal terms involving the general 
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T) nm (u>o) with 9 n ^ 6 m . However, if the transformation matrix is such that the 
sums over 9 n , 9 m are dominated by the contribution coming from a single (real) 
diagonal element D„„ with transverse mode index 9 n , then the other terms can 
be ignored, leading to the following simple result for the decay rate: 



= i^rf, (146) 



where 



rf = 4np(Lu )Y,\9n(^8n)\ 2 (147) 



is the decay rate expected for a normal cavity situation. Thus we see for this 
special case where a single NHM dominates, the decay rate has been enhanced 
by the Petermann factor K n . 

4 Conclusion 

The present paper is a fully quantum treatment of the field for unstable optical 
systems. A standard canonical quantization proceedure is used, based on ex- 
panding the vector potential in the unstable cavity region and in the external 
region via non-Hermitean (Fox-Li) modes and their Hermitean adjoint modes, 
all defined via the optical system. Three dimensional systems are treated, us- 
ing the paraxial fl22fl and monochromaticity approximations |19| for the cavity 
non-Hermitean modes. Both right and left travelling modes are included. The 
results are similar to JliJ, |2l[] but differ in detail. The field is equivalent to a 
set of quantum harmonic oscillators (QHO), associated now with non-Hermitean 
modes rather than true modes, and thus confirming the validity of the photon 
model for the case of unstable optical systems. The annihilation, creation oper- 
ator pairs for each QHO are not Hermitean adjoints. A doubling of the number 
of annihilation, creation operator pairs occurs compared to |19|| , |2lf| , showing 
that the total number of true mode QHO's equals the total number of QHO's 
associated either with the cavity or the external region non-Hcrmitcan modes. 
The final form of our quantum Hamiltonian for the EM field is the sum of non- 
commuting cavity and external field contributions. To a good approximation, 
both the cavity and external field Hamiltonians can be expressed as a sum of 
independent QHO Hamiltonians for each non-Hermitean mode, but the exter- 
nal field Hamiltonian also includes a coupling term responsible for external non- 
Hermitean mode photon exchange processes. The two independent QHO Hamil- 
tonians are alternative choices for an unperturbed Hamiltonian. Left and right 
eigenstates for each of these unperturbed Hamiltonians are obtained, and the 
energy is given by the usual QHO result. Certain cavity non-Hermitean mode 
annihilation, creation operators do not commute with others for the external 
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region non-Hermitean modes, leading to cavity energy gain and loss processes. 
A simple description of the cavity-external region light coupling in terms of sur- 
face integrals involving products of cavity and external non-Hermitean mode 
functions is found. 

Atom-field interactions are treated in the electric dipole and rotating wave 
approximations. Atomic transitions leading to cavity non-Hermitean mode pho- 
ton absorption are associated with a different coupling constant to that for 
atomic transitions leading to photon emission. This feature is directly related 
to the treatment being based on non-Hermitean mode functions and leads to 
enhanced emission rates. Using the essential states approach, the spontaneous 
decay of a two level atom located in the cavity is treated. The external region 
term in the field Hamiltonian is neglected assuming that atomic decay into the 
cavity is much faster than cavity decay. Coupled equations for the amplitudes of 
atom-field states involving the excited atom and no photons or the ground state 
atom with one photon in a cavity non-Hermitean mode are obtained. Marko- 
vian decay occurs under certain conditions, the decay rate being enhanced by 
the Petermann factor p2| in special cases when a single NHM dominates. 
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6 Appendix A: Commutation rules for cavity 
and external region operators 

6.1 Commutation rules for D and B field operators 

To derive the commutation rules between cavity NHM generalised coordinates 
or momenta and external region NHM generalised momenta or coordinates we 
use the general rules for the magnetic field and electric displacement operators: 

[B a (R),D /? (R / )] = ^£ Q/37 9/9A 7 (5(R-R / ) (148) 

[D Q (R),B /3 (R / )] = -ihea/^d/dX^S^R-B/), (149) 

where R is in the cavity region and R' is in the external region. In the cavity 
region we will use the lowest order terms in the small parameter / for both 
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B(R) and D(R) namely: 

B(R) = ^(ifc„Q„zxU„(R)- i fc nJ R] l zxV;(R)) (150) 

D(R) = -^(5„U„(R)+P r tv;(R)), (151) 

n 

whilst in the external region we will write: 

D(R/) = -^(5kU^(R / )+^V|,(R/)) (152) 

K 

B(R/) = ^(feVxU K (R/)+4VxV^(R/)). (153) 

K 

For simplicity the cavity will be taken to be the region between the planes z = zq 
and z = zi, + e and the external region to be the region between z = zj, — e and 
z = oo. The quantity e is a small distance to allow for the cavity and external 
regions to overlap slightly, and will be taken to zero afterwards. The plane 
z = zjj is the boundary between the cavity and external regions, which are 
shown in figure 2. 

6.2 Cases of cavity NHM coordinates and external NHM 
momenta 

The commutation rules between cavity NHM coordinates and external re gion 



mom enta are found by substituting for B(R) and D(R/) from equations (150, 



152 ) into equation (148), which gives: 



ih J2 e «/37 d/dX-i <KR - R 7 ) = - J2J2( lkm eazjUmCR) u[ (R 7 ) [Q m , S L ] 

1 7 m (7) L 

+ ik m e az -yU m (R) Vf *(R/) [Q m , PI] 

-ik m e azi V^{TL) U[ (R/) [Rl, S L ] 

- ik m e«WC(R) Vf *(R/) [Rl, P[] }. 

(154) 

On the right side the sum over 7 only involves x, y and that over m only those 
with & m = •y. For such terms the components of U m and VJ^ are U m and V* n . 



6.2.1 Case of [Q„, P i 



K\ 



To obtain the commutation rule for [Q n , P K ] we multiply each side of equation 
(154) by V*{R) f/^(R/), integrate R, R' over the cavity and external regions 



respectively, then sum over [3. Using the biorthogonality results in equations 
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LOL GO) and neglecting integrals where the complex conjugation occurs either 
zero or two times we find that: 

ih V£a/3 7 / d 3 R [ d 3 R/ V*(R) U K {R')d/dX^5{R- R f ) = - V] ik n e azj [Q n , P} { ]. 

M Jc Je 7=«» 

(155) 

There are only two non-zero possibilities a = x or a = y. These give the two 
equations: 

ik n {Qn, Pr]- = ih [ d 3 R f d 3 R/y„*(R) U^CR^d/dzSCR-R,/) 
Jc Je 

-ih f d 3 R I d 3 R/V*(R) Ul c {R')d/dy8{R-R / ) 
Jc Je 

(156) 

with a n = y and 

- ik n [Q n , P^]. = ih f d 3 R I d 3 R^V*(R) U z K {R / )d/dx5{R-R / ) 

Jc Je 

-ih f d 3 R f d 3 R/v;*(R) UK(R / )d/dzS(R-R / ) 
Jc Je 

(157) 

with & n = x. 

Considering the double integral on the right side of equation (157) involving 
d/dx, this is of the form: 



f dydz [ d 3 R/ 5{y-yl) ^{z-z 1 ) U^R') 
Jc Je 

dxV*(R)d/dxS(x-x / ) 



I x = 

IE 



(dydz f d 3 Rl 8(y-yl)8{z-zl)U z K {Rl 
Jc Je 



■ ■•■•- ■"■■■< ■■< • "•• - - ' • i\ • ■•■' 
ic Je 

dx{d/dx[S(x-x/)V*(R)}~S(x- x /) d/dx [V*(R)]}, 
c 

(158) 

using partial integration. For the first term carrying out the x integration leads 
to a result dependent on V*(R) for x — — oo and x = +oo. As the NHM 
functions eventually go to zero for large x, this contribution can be ignored. For 
the remaining term we now have: 



I x = 



/ d 3 R f d 3 R/<5(R-R / ) 17£(R/) d/dx[V*(R)} 
Jc Je 

d 2 s f dVtf(s-V) U z K {s',z b ) d/dx[V*(s,z b )] 

< Je 

z b +i roc 

dz / dz' Siz-z 1 ), (159) 

ZO J Z b —£ 
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where in view of the <5(R — R' ) factor, the ordinary functions can be expressed 
in terms of their values on the boundary, where z = Zb- The integrals over z and 
zl give a result proportional to the small quantity e, and thus the integral I x 
is zero. Similar considerations show that the term involving d/dy on the right 



side of equation fll56| ) is also zero. The two equations fll5q , 157) are now both 
the same as: 

[Q n ,Pp\. = (fi/Av,) / d 3 K f d 3 R/v;(R)- Ujr(R/)0/0*tf(R-B/) 

JC J E 

(160) 

since for the case where a n = x the scalar product picks out the component 
!7jf(R') and for the case a n — y the scalar product picks out the component 
0*(R/). 



Considering the double integral on the right side of equation ( 160 ) involving 
d/dz, this is of the form: 



h = 

ic 



f d 3 R f d 3 R/v;(R) • Ujf(R/) d/dz [<J(R- R 7 )] 
= / dxdy f d 3 R/ <5(x-x / )'5(2/-2/ / )Uk(R / ) 

• / d«{0/a?[<5(*-«/)V*(R)] -*(«-«/) <9/az[V;(R)]} 
ic 

= / d 2 s /" dV <5(s - s 7 ) / dz 7 U^s 7 ,^) 

JC ifi -izt-e 

• /" dz { <9/<9z [<5(*-a/)V* (s,«)] - (5(z-z/) <9/<9z [V*(s,z)]} 

izo 

(161) 

using partial integration and then specifying the limits of integration over z and 
z< , The integration over z in the first term is carried out and unlike the similar 
integrations over x or y that applied to the I x or I y cases we find that a non-zero 
contribution results from the z = Z]> + e limit. The second term is similar to the 
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other term in the I x or I v cases and overall we now have: 



I, = 



f d 2 s [ dV S(s - s/) f dz 1 \3 K (*I ^z 1 ) 

JC JE Jz b -e 

■{S(z b + e-z^V^fazb + e)} 

- [ d 2 s f dV <5(s - s 1 ) t dz 1 Uk^ 1 ,* 1 ) 

JC JE Jz h -e 

■I' dz{5(z-zl) d/dz[V* n {s,z)}} 

J z 

[ d 2 s f dV<5(s-s/)V*(s,,? 6 ) 

JC JE 

/>oc 

■Uk(sU) / dz 1 5(z b + e-zl) 

J Zb — e 

- / d 2 s / dV <5(s - s /) Ujf(s/, z b ) • 5/dz [V* (s,z)] Z6 



A' 
oo 



«fa/ 



Zb+e 



dz 5(z~z') 



(162) 



In equation (162) the first term is non-zero, since the z' integral equals one and 
we are left with a surface integral over the cavity-external region boundary when 
the 5(s — s/) factor is integrated out. The second term in equation ( |162| ) is zero, 
the integrals over z and z' giving a result proportional to the small quantity e, 
and as in the I x case the overall integral is zero. Hence we have: 



I z = I d 2 sV* n (s,z b )-U K (s 7 z b ) 
Js 



and thus 



[Q„, PT.]. = (h/k n ) / d 2 s V:(s,z b ) • V K (s,z b ). 

JS 



(163) 



(164) 



6.2.2 Case of [R n , & K ] 

To obtain the commutation rule for [R n , S K ], we first obtain the result for 
[R'l, Sk] by multiplying each side of equation (154) by f/„(R) V£*(R/), inte- 
grating R, R' over the cavity and external regions respectively, then summing 
over /3. The treatment then involves the same steps as in the case of [Q n , Pr] 
except for the replacement of V* by U n and U^ by V£* toget her w ith the pres- 
ence of a negative sign resulting from the original equation (154) itself. This 
leads to: 



[RlS K ]=-(h/k n ) [ d 2 S U n ( S ,z b )-V* K (s,z b ), 
Js 



(165) 
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so after taking the adjoint of both sides of this result we have finally: 

[R n ,S\ c } = (h/k n ) /d 2 sU,*(s,z h ).V K (s,2 6 ). (166) 

Js 

6.2.3 Cases of [Q n , S K ] and [R n , Pj { ] 

To obtain the commutation rules for [Q n , S K ] and [R n , P K ] we use equations 
(Ell) to write Q n as a linear combination of the R m and R n as a linear com- 
bination of the Q m . This involves the transformation matrices D nm and C nm 
respectively. We then use the previously derived results in equations (166|, 164) 



for the commutators [R m , S K ] and [Q m , P K ]- The (h / k m ) terms in the sums 
over m can be replaced by {h / k n ) since the transformation matrices are zero 
unless k n = k m . We then use equations (jlj) and the Hermiteancy properties of 
D nm and C nm to perform the sum over m, and this leads to the NHM functions 
U^j and V^j being replaced by V* and U* respectively. This gives the results: 

[Qn,& K ] = (&/*») /d 2 sV;(s,z 6 )-V K (s,z 6 ) (167) 

Js 

[Rn,P J K ] = (h/k n ) [ d 2 sXJ* n (s,z b )-XJ K (s,z b ). (168) 

Js 

6.2.4 Cases of [Q n , P K ] , [R n , S K ],[Qn, S K ], [Rn, Pk] and their ad- 
joints 

These commutators involve either two or zero adjoints and we conclude they 



may be taken as equal to zero. Starting from equation ( 154 ) we could multiply 
each side by V^(R) V£*(R/), integrate R, R' over the cavity and external re- 
gions respectively, then sum over (3 and follow similar steps as for [Q n , P K ] to 
obtain a surface integral expression for [Q n , Sk]- The surface integral would 
involve the scalar product V*(s,Zb) • V* K (s,Zb) with both factors conjugated. 
The cavity-external boundary is never really planar, so the exponential factors 
associated with V*(s,Zfc) • V^ (s,Zb) would average out to zero. On this basis we 
neglect commutators such as [Q n , Sk], [Rn, Sk] , [Qn, Sk], [Rn, Pk] and their 



adjoints. Alternatively, using equations (18) for the cavity NHM generalised 
coordinates and those analogous to equation ( |4q ) for the external NHM gener- 
alised momenta, we see that such commutators would be linear combinations of 
true mode commutators [qk,Pi], [<lk^ ,p\] that are all zero. 

6.3 Cases of cavity NHM momenta and external NHM 
coordinates 

The commutation rules between cavity NHM momenta and external regio n co- 



ordi nates are found by s ubstituting for D(R) and B(R/) from equations ( 151 



153|) into equation (149), which gives: 
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7 m(a) L 

+U m (R) {V x V L (R/)}' 3 * [S m , R{] 
+C(R){VxU L (R/)f[Pi,Q L ] 

+y:(R){Vxv L (R/)f[AU[]}- 

(169) 
On the right side the sum over m only involves terms with a m = a. 

6.3.1 Case of [S n , R ] K ] 

To obtain the commutation rule for [S n , R K ] we multiply each side of equation 



( [J69l ) by V*(R) {V x Ux(R / )} /3 , integrate R, R 7 over the cavity and external 
regions respectively, then sum over (3. Using the biorthogonality results in 
equation (nil) and the results in equations (B9, Eoh for the magnetic energy term 
and (as usual) neglecting integrals where the complex conjugation occurs either 
zero or two times we find that: 

k 2 K [S n , R f K ] = ih V £ Q/37 / d 3 R / d 3 R/ V:(R) {V x U* (B/)} /3 

7 JC J E 

. 8/dX-, d(R - R/) 

(170) 

where a„ = a. Since a„ = x or a n = y are the only two possibilities we get 
the two equations: 

fc| [S n , #Jj.] = ih [ d 3 R [ d 3 R f V*(R) {V x U K (R / )} 1 '5/az(5(R- R 7 ) 

JC J E 

-ih [ d 3 R [ d 3 R/v*(R){V xU k (R/)} z d/dySiR-R/) 

JC J E 

(171) 
for a. n = x and 

fc| [S n , AJj-] = - ih [ d 3 R [ d 3 R/ V*(R) {V x U x (R/)} x a/az(5(R-R / ) 

+ ih J d 3 R J d 3 R/v*{R) {V x U K (R/)} Z d/dx 5(R- R/) 

JC J E 

(172) 

for a n = y. 

As previously, the integrals involving d/dx or d/dy give zero, so only the 
integrals involving d/dz remain. Following the same treatment for the latter 
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results in expressions involving surface integrals: 

k 2 K [S n , R K ] = ih [ d 2 s V:(s,z b ) {V x V K (s, z b )}y (173) 

Js 



for OL n — x and 



k 2 K [S n , R K ] =-ih [ d 2 s v:(s,z b ){V x V K {s,z b )} x (174) 

Js 

or cx n = y. These results may be summarised in a single formula: 

[S n , R K ] = (ih/k 2 K ) f d 2 s z ■ V* n (s,z b ) x {V x U K (s,z b )}. (175) 

Js 

6.3.2 Case of [P n , & K ] 

To obtain the commutation rule for [P n , Q K ], we first obtain the result for 
[P,\, Qk] by multiplying each side of equation ( fL69| ) by U n (R) {V x V K (R/)}^*, 
integrating R, R' over the cavity and external regions respectively, then sum- 
ming over p. The treatment then involves the same steps as in the case of 
[S n , PJ K ] except for the replacement of V* by U n and {V x U^ (R/)}' 3 by 
{V x V K (R / )} /3 *- This leads to: 

[PI Q K ] = {ih/lsfc) f d 2 sz-XJ n (s,z b ) x {V x V^(s,z fc )}, (176) 

Js 

and taking the adjoint of both sides leads to the result: 

[P ni $ K ] = (ih/k 2 K ) f d 2 s z ■ U* n (s,z b ) x {V x V K (s,z b )}. (177) 

Js 

6.3.3 Cases of [P n , R K ] and [S n , Q K ] 



To obtain the commutation rules for [P n , R K ] and [S n , Q K ] we use equations 
(E9J) to write P n as a linear combination of the S m and S n as a linear combination 

K ] and [R n , PI] 



of the P m . The same method used for determining [Q n , S K ] and [Rn, P^] is 



then followed and we find that: 



[P n ,R K ] = (ih/kD frf 2 sz-U;M x{Vx\J K (s,z b )} (178) 
Js 

[S n , Q K ] = (ih/k 2 K ) f d 2 s z ■ V n (s,z b ) x {V x V K (s, z b )}. (179) 
Js 

6.3.4 Cases of [P n , Qk] , [S„, Rk],[Pu, Rk], [S n , Qk] and their ad- 
joints 

These commutators involve either two or zero adjoints and we conclude they 
may be taken as equal to zero. The details are similar to those for [Q„, Pk] , 
[Rn, Sk] , [Qn, Sk], [Rn, Pjc\ and their adjoints. 
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7 Figure captions 

Figure 1. Unstable optical resonator system consisting of two confocal mirrors. 
The dashed line indicates the boundary between the cavity and the external 
regions. The left hand mirror is assumed to be large, so that light does not 
travel past it and is confined to the cavity and external regions shown. 
Figure 2. Ranges of z integration for cavity and external regions used to deter- 
mine cavity NHM/external NHM commutation rules. 
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